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Abstract 

We discuss various definitions of decoherence and how it can be measured. We compare 
and contrast decoherence in quantum systems with an infinite number of eigenstates (such 
as the free particle and the oscihator) and spin systems. In the former case, we point out the 
essential difference between assuming "entanglement at all times" and entanglement with the reser- 
voir occuring at some initial time. We also discuss optimum calculational techniques in both arenas. 
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I. INTRODUCTION 



There is presently intense interest in the fundamentals of quantum theory and appli- 
cations. In particular, the superposition principle (and related work on Schrodinger cats), 
entanglement, and the quantum-classical interface are at the cutting-edge of topical research, 
especially in respect to their relevance to quantum computing, quantum information pro- 
cessing, quantum teleportation and quantum encryption. Since superposition states are very 
sensitive to decoherence, reservoir theory has attracted much recent interest. 

Decoherence is the physical process which is responsible for the emergence of the classical 
world from the quantum world . In order to understand the phenomenon 

in depth, much attention has been devoted to the case of a free particle interacting with 
a reservoir. Most of these investigations have assumed that entanglement with the reser- 
voir occurs at some initial time, leading to the conclusion that a dissipative environment 
is necessary to achieve decoherence. However, we have shown that in the case of "entan- 
glement with the reservoir at all times" , decoherence can occur simply due to temperature 



effects without requiring dissipation 



10[. In addition, we pointed out that there are 



many definitions of decoherence and we argued that the preferred definition involves 
probabilities in coordinate space since they can be measured (in contrast to the off-diagonal 
components of the density matrix, as considered often in the literature). Moreover, the fact 
that the free particle wave-packet spreads even in the absence of a reservoir necessitates a 
careful consideration of how decoherence should be defined. 

By contrast, an oscillator wave-packet does not spread in the absence of a reservoir so 



that its investigation was clearly worthy of investigation 



12| . This brought to light another 



interesting feature viz. the fact that revivals of coherence can occur. A new dimension was 



added to the problem by consideration of an external field 



15 



3. 



13 



1J|, of interest because of 



the recent experiments of the Wineland group 

Thus, in Sees. II and III, we present a detailed discussion of decoherence for both a free 
particle and an oscillator, respectively, in an arbitrary reservoir, including remarks on what 
we consider to be the optimal tools for such calculations. Sec. IV is devoted to a discussion 
of decoherence in spin systems and the significant differences that arise compared to systems 
with an infinite number of eigenstates. In Sec. V, we summarize and discuss our conclusions. 
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II. DECOHERENCE FOR A FREE PARTICLE 



Decoherence refers to the destruction of a quantum inferference pattern and is relevant 
to the many experiments that depend on achieving and maintaining entangled states. Ex- 
amples of such efforts are in the areas of quantum teleportation, quantum information and 
computation, entangled states, Schrodinger cats, and the quantum-classical interface. Much 
of the discussion of decoherence has been in terms of a particle moving in one dimension 
that is placed in an initial superposition state (Schrodinger "cat" state) corresponding to 
two widely separated Gaussian wave packets, each with variance and separated by a 
distance d. We consider the particle to be coupled to an arbitrary reservoir such that in the 
distant past the complete system is in thermal equilibrium at temperature T. For such a 
state the probability distribution at time t can be shown to be of the form 



= Pi + P-2 + 2P,case(t) (2.1) 
where is the probabihty distribution for a single wavepacket, given by 

Po(a;,t) = , ^ exp|--47-|. (2.2) 
Here and in (j^.lj) w'^{t) is the variance of a single wavepacket, which in general is given by 

w\t)=a'-''^^^^^0^ + s{t), (2.3) 
where cr^ is the initial variance, [x(0),x(t)] is the commutator, and 

s(t) = {{x{t) - xm") , (2.4) 

is the mean square displacement. The temperature dependence enters only in s{t). Also, 
a(t), which can be defined as the ratio of the factor multiplying the cosine in the interference 
term to twice the geometric mean of the first two terms |8] is given by the following exact 
general formula 



3 



a{t) = exp 



(2.5) 



8aV{t) J ■ 

The question now arises as how to define decoherence. It is often described as the disap- 
pearance of the interference term with time. However, this is not strictly correct since the 
integrated probabihty of each of the three terms in ()2.H1 is constant in time. In fact, if we 
define the common normalization factor in ()2.1|) as 



N 



2 (l + 



-1/2 



(2.6) 



we obtain 



r dxP{x, t) = AT^ i 1 + 1 + 2 exp 

■J—OO 



■2 



^2iV^|l + exp^ 8aVJ 
= 1- (2.7) 

It is clear that the integrated probability of each of the 3 terms is constant in time. Also, the 
integrated probability of the interference term is smaller than the integrated probability of 
either of the direct terms by a factor exp (— ^) which, for macroscopic d and microscopic 
a is very small. Thus, for example, when (rf/cr) = 5 we see that exp {^—-^^ ~ 4.4 x 10~^ 
which is a very small number compared to unity. 

In fact, the effect of decoherence is to cause the interference wave packet to spread in time, 
with a concomitant decrease in its amplitude, consistent with the fact that its integrated 
probability is constant in time. Moreover, we note the ubiquituous cos term which is a 
feature of the interference contribution. Thus, a measure of decoherence is to focus on the 
interference term by calculating and measuring the time dependence of Pj. However, if 
decoherence is taken to mean the effects of a reservoir (not necessarily dissipative) then this 
definition has a problem since it given a non-zero value for the rate of decoherence in the 
absence of a reservoir. The problem stems from the fact that the wave pac ket for a free 
particle spreads naturally in time, a standard result in quantum mechanics jl2|. In fact, this 
result immediately follows from ()2.3|) since, for a free particle without dissipation and with 
zero temperature, s{t) = and [a;(0), a;(t)] = iht/m, so that 
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2+2 



(2.8) 



As a result, all of the three wave packets in 1)2.11) spread in time. Thus, in order to correct 
for this dynamical spreading, we were led to define 



a(t) 



Pi 



(2.9) 



as a measure of decoherence. It has the virtue of reducing to unity in the absence of a 
reservoir [lol |. 

In the case of Ohmic dissipation, high temperature {kT » h'y), small times {'-ft << 1) 
and d » Xth, cr, we obtained 0,1^ s{t) = {kT /m)t^ and w'^{t) a^. Hence ()2.5|) reduces to 



a(t)^exp|-(l) }, 



(2.10) 



where 



(2.11) 



'kT/m d 

exhibiting "decoherence without dissipation". Here 7 is the Ohmic decay rate, and Xth 
h/^/mkT is the deBroglie wavelength. Hence 



A?,.^(5.2X10-V,„)=(1^)(5^). (2,12) 

For example, the choice m = Igm, d = 1cm and T = 3007^, leads to the very large ratio 
d/Xth = 2 X 10^°. In addition, we note that 

^ = (2 13) 

h'y 7(10iis-i)' n 

On the other hand, if we consider a low temperature environment {kT << hj) [l^, then 
using the quantum Langevin equation j3| to evaluate w'^{t), we obtained 



a{t) 



exp <; \ - 



, (t 3 
log — + 7ij - - 
m 2 



|, r«t«(C/m)-i 



(2.14) 



with 



_ ma 6TT 



(2.15) 
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and ( = 7m and where 'Je = 0.577215665 is Euler's constant. Thus, in the low-temperature 
regime, decoherence requires a dissipative environment. 

We emphasize again that all of the above results correspond to the case of the free 
particle being entangled with the environment at all times. Initially (or in the distant 
past), the complete system is in equilibrium at temperature T |^]. At time t = 0, say, a 
first measurement is carried out which prepares the system of two widely separated wave 
packets. Then, after a time t, a second measurement is carried out which probes the system. 

Most other investigations have assumed that the quantum system is initially at tem- 
perature zero and uncoupled from the reservoir, which is at temperature T. When the 
quantum system is brought into contact with the reservoir at t = 0, it takes a time ~ 7"^ 
in order for thermal equilibrium to be achieved, which is generally much larger than the 
decoherence decay time so that decoherence occurs before thermal equilibrium is achieved. 
Such a calculation generally involves the use of density matrix equations. Thus, in order to 
make a careful comparison with our results corresponding to "entanglement at all times" 
we solved the exact HPZ master equation |3] and obtained results for an arbitrary reservoir 
and arbitrary temperatures. In the particular case of high temperature, we obtained 

If we suppose that the slit width is negligibly small, we find a(t) = exp{— t/r^i} where 
Trf = ^H^- This is similar to the decoherence time that often appears in the literature. 
But, as we have mentioned above, this result corresponds to a particle in an initial state 
that is effectively at temperature zero, which is suddenly coupled to a heat bath at high 
temperature. The result is therefore unphysical in the sense that the initial states does not 
correspond to that envisioned when we speak of a system at temperature T. Thus, with this 
scenario, we get no decoherence, [i.e. a{t) = 1] when 7 = 0. On the other hand, we found 
that we could repair this unphysical difficulty by choosing the initial temperature to be the 
same as the reservoir, in which case the result given in ()2.1()j) is again obtained. 

We conclude this section by commenting on the techniques used. The "entanglement at all 
times" calculation utihzed quantum probability distributions (which are related to Wigner 
distributions which are probabilities in quantum phase space 



21 



22| |) in conjunction 



with results obtained by use of the stationery solution to the generalized quantum Langevin 
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equation [igj ]. 

For the case where the oscillator and reservoir are initially decoupled, we formulated the 
problem in terms of the Langevin equation for the initial value problem. In fact, a master 
equation was not required but, in order to make contact with other investigators, we actually 
derived the HPZ exact master equation from the Langevin equation approach. Moreoever, 
the Langevin approach enabled us to obtain explicit results for the time-dependent coef- 
ficients in the HPZ equation. Our strategy was based on use of the Wigner distribution 
eventually leading to explicit and very general results for the Wigner distribution for the 
reduced system of an oscillator in an arbitrary state and in an arbitrary heat bath. These 
results were then used to obtain coordinate probabilities. 

The Wigner distribution is, of course, the Fourrier transform of the density matrix |2ol |. 
so that results obtained by use of the former are equivalent to results obtained by use of the 
latter. However, Wigner functions are much easier to use in practice because, as distinct 
from density matrices, they are not operators and they are always real. They describe, 
in essence, quantum phase space. However, just as we object to using the density matrix 
as a quantitative way to describe decoherence, we also would have the same objections to 
the use of Wigner distributions for that purpose ^| but, instead, we regard them as a 
wonderful tool to calculate probability distributions, which are measurable quantities. A 
detailed review of the use of Wigner distributions in this context appears in j2l|. 

III. DECOHERENCE FOR AN OSCILLATOR 

A free oscillator wave packet does not spread in time (i.e. its width does not depend on 
time) but it oscillates back and forth with a frequency cj, where uj is the oscillator frequency 



22| . In other words, its shape does not change but the peak of the wave packet has the time 
dependence xq cosut. Thus, in contrast to the free particle, there is no dynamical spreading 
for an oscillator wave packet; it is a coherent state. 

In the case of a Schrodinger cat superposition of two coherent states separated at a 
distance d and in thermal equilibrium in a non-dissipative reservoir, we obtained a result for 
P{x, t) which is similar in structure to the result given in ()2.Hl . the difference being that the 
interference term, in common with the other terms, is coherent in the sense that the wave 
packet shape does not change but, instead, it oscillates in time and persists for all time. 
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Explicitly, the corresponding result for a{t) is 

, , I muj<P cos^ Lot . . . 

a[t) = exp < } . (3.1) 

[ 2/isinli(f^) ' 

It is interesting to note that for small times {ut << 1) after the times for which the 
attenuation factor has its maximum value of unity (cos tut = 0), and if we recall that the 
initial width a of the oscillator wave packet is given by = {h/2muj), then for uj 
we find that a{t) reduces to the result given in ()2.1()|1 and ()2.1H1 for the free particle case. 
However, in contrast to the latter case, for later times, there is a revival of a{t) toward 
its maximum value. We believe that this revivial of coherence in a non-dissipative thermal 
reservoir is a common feature, the only exception being the free particle. For the case of 
7 7^ 0, we refer to 

Finally, we remark that a non-random external field f(t) does not cause decoherence, in 

nn 

contrast to the case of a random / (t) [Ij, ll^ , a result of relevance for the analysis of the 



experimental results obtained by the Wineland group 



IV. DECOHERENCE IN SPIN SYSTEMS 
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Whereas the study of decoherence for a free particle and an oscillator is very relevant for 
the study of the classical-quantum correspondence, it is now generally believed that spin 
systems (where the emphasis is on the spin and associated magnetic moment of quantum 
particles) are more relevant to quantum computing, teleportation and information process- 
ing. There already exists extensive studies of such systems 24, 3,13 confine ourselves 



to some general remarks. 

The main advantage of spin qubits are that they interact weakly with their environment. 
In general, as distinct from the study of free particle and oscillator systems interacting with a 
reservoir, we do not have to worry to the same extent with the question of when entanglement 
with the environment occurs since the system is generally controlled by external fields. Thus, 
density matrix techniques will invariably constitute the tool of choice. This brings up the 
question as to whether or not there is an optimum way to solve these operator equations. We 
already saw in Sees. II and III that Wigner distribution functions proved to be a wonderful 
calculational tool for the study of system with an infinite spectrum of states. However, 
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for spin systems, there is also what we consider is a very useful and physically appealing 
tool viz. the use of the spin polarization vector. This has found application for particles 
of arbitrary spin in a magnetic field but not interacting with a reservoir 27j but, strangely 
enough, it has not usually been adopted for the more challenging situations where a reservoir 
is present. Ford and me actually used this technique in order to calculate the spectrum of 
resonance fluorescence for a driven two- level atom instead of using the corresponding 
but more unwieldy master equation. Thus, in the Appendix, we use the spin polarization 
vector to study an even simpler system. This is the well-known two-level system (describing 
either a spin | system or a two- level atom) interacting with a reservoir and we feel that, our 
calculation demonstrates the simplicity and physically appealing nature of using the spin 
polarization vector. 



V. DISCUSSION 



We have pointed out that, whereas it is easy to describe qualitatively what is meant 
by decoherence, difference of opinions may arise as how best to define its signature in a 
quantitative manner. Different results may ensue depending on 

(a) whether one is dealing with coordinate, momentum or quantum phase space probabilities 
(and we favor the former since they are measurable). 

(b) how natural dynamical spreading of a wave packet in the case of a free particle is 
separated from spreading due to environmental effects. 

(c) whether one considers entanglement with a reservoir as existing for all times or simply 
taking place at some initial time. 

(d) whether external forces come into play. 

(e) the system under consideration since quantum systems with an infinite number of 
eigenstates (such as the free particle and the oscillator) generally exhibit decoherence decay 
times much smaller than relaxation times whereas, for spin systems, decoherence decay 
times are often comparable to relaxation times (which is the reason they are preferred for 
applications). 

In addition, we pointed out the great calculational advantages of using either Wigner 
functions (for the free particle and the oscillator) or spin polarization vectors (for spin 
systems) . 
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APPENDIX A 



Solution of the master equation for a two-level system interacting with a 
reservoir, using the spin poleirization vector. 

The master equation for the reduced density matrix p may be written in the form 

1 

= 2^^^^^ (2(7_p(7+ - (7+(T_p - p(7+(7_) 

1 - 

+ 2^ (2cx+pcT_ - a-a+p - pa-a+) , (Al) 
where n is the thermal average boson number 

n = [exp {hcu /kT} - 1] ( A2) 

Also 

o-± = ^ {c^x ± iCTy) , (A3) 

where the (jj are the usual Pauli spin matrices. 

Now p can be expanded in terms of the complete set (/, (7^) with real coefficients: 



p^^(l + p.a), (A4) 

— * 

where / is the identity matrix and P is the polarization vector. This leads, in the usual 



manner to 

P={a)=Triap), (A5) 
i.e. the polarization vector is the ensemble average of the spin vector. Thus 

Next, we note that 
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Tr {a {a^a^p + pcr+cr- — 2(T_pcr_|_)} 
= (cr(T_|_cr_ + cr_|_o"_a — 2(T_|_(Tcr_) 

(1 + a,) + ^ (1 + a,) a + (1 + a, 
((J + z + i + (T^z) 
P + z ■ Pz + 2z. 



and 



Tr {(? (cr_(T_|_p + p(T_(T+ — 2(T_|_p(j_)} 

= {a — 2z + GzZ) 
= P + z ■ Pz ~2z. 

Hence, using (EH), and (jMl) to (IXsjl . 



It 



At equilibrium (where P = Pq) 



so that 



Hence, the right-side of ()A9|1 is 



-ln[P + PJ-2z) 
-l{2n+l) ^p + PJ^-^z. 



dt 



dt 



zero, which imphes 

P + P,5: 



2n + 1 
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^ — * 

It is thus clear that, in thermal equilibrium at temperature T, P is along z i.e. 



P = P,z = PqZ. (A13) 

Thus, from ()A12|) . we obtain 

Po = — = - tanh { ^\ . (AM) 

For T = 0(n = 0), we see that Pq = —1 corresponding to all the spins being down. On the 
other hand, for T —>■ oo{n oo), we obtain Pq —>■ 0, corresponding to the equal number of 
up and down spins. Hence, ()A9|) becomes 

^ = 7^(P + PJ)-7Z. (A15) 



dt 2P 



Hence 



dt Po " ^ 

Po-Pz 



(A16) 



where 



1.-^ = 0(2« + I)=ocoth(|^). (A17) 

Thus, for T — > 0, Ti ^ 7^^ whereas for T — >^ 00, Ti — > 0. The solution of ()A16|) is 

n = n{l-exp(-^)}. (A18) 
In addition, ()A15|) may be written as 

dP 1 

^P + P,z) - -fz. (A19) 



dt 2Ti 

Thus, defining Pi to be the component of P perpendicular to z, we obtain 



(A20) 



dt 2Ti 
whose solution is 
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P± = P±(0)exp(-^), (A21) 



where 



T2 = 2Ti. (A22) 
Since the magnetic moment M of a nuclear ion with spin S = f ci is given by 



a e ^ eh 

M = Qn o = a 

2mc 2mc 2 

gn ^ 

where m is the nuclear mass, e is its charge, Qn is the nuclear ion g factor and 



2mc 

is the nuclear Bohr magneton, it follows that 



2 N / 2 

Hence, using (IA18|) and ()A21|) . we obtain 



and 



(A23) 



/io = (A24) 



m) = -f^o^ (<?) = -^o%A (A25) 



(M,) = M„|l-exp(-^)}, (A26) 



(Mx) = M^(0)exp(-^). (A27) 

We note that Ti is the so-called spin-lattice (longitudinal) relaxation time which describes 
the approach to thermal equilibrium whereas T2 is the dephasing (transverse) relaxation 
time, describing a system of spin | particles in a magnetic field. 
It is also of interest to re-write ()A4|1 in the form 

p = ^ {1 + (a+P_ + a_P+) + a,Pj , (A28) 

where 
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P± = (P, ± iPy) 

SO that, as is also obvious from ()A19|1 to ()A2H) . 



(A29) 



and 



Moreover, if we write 



and use the fact that 



dP+ P4 



dt 2Ti 



P± = P±(0)exp(-^ 



it readily follows from (jA28|) that 




1 





(A30) 



(A31) 



(A32) 



(A33) 



(A34) 



1-P. 



(A35) 



and 



(A36) 



(A37) 



Thus, using (1A18|) and ()A3Hl . it is clear that we obtain the familiar result that the rate of 
decay of the diagonal elements of the density matrix (^r^) is twice that of the rate of decay 
of the off-diagonal elements (^2^^) • However, we suggest that it is physically more appealing 
to say that we conclude [from ()A18|) and either ()A21|) or ()A31|) ] that the expectation value 
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of the spin component in the z direction (direction of the magnetic field) decays faster than 
the corresponding values perpindicular to the z direction. 
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